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We apply the technique of parameter splitting to existing cosmological data sets, to check for a generic failure 
of dark energy models. Given a dark energy parameter, such as the energy density S7a or equation of state 
w, we split it into two meta-parameters with one controlling geometrical distances, and the other controlling 
the growth of structure. Observational data spanning Type la supernovae, the cosmic microwave background 
(CMB), galaxy clustering, and weak gravitational lensing statistics are fit without requiring the two meta- 
parameters to be equal. This technique checks for inconsistency between different data sets, as well as for 
internal inconsistency within any one data set (e.g., CMB or lensing statistics) that is sensitive to both geometry 
and growth. We find that the cosmological constant model is consistent with current data. Theories of modified 
gravity generally predict a relation between growth and geometry that is different from that of general relativity. 
Parameter splitting can be viewed as a crude way to parametrize the space of such theories. Our analysis of cur- 
rent data already appears to put sharp limits on these theories: assuming a flat universe, current data constrain 
the difference AOv = f2A(geom)-f2 A (grow) to be -0.0044^:m57^:oii9 (68% and 95% C.L. respectively); al- 
lowing the equation of state w to vary, the difference Aw = w(geom) - w(grow) is constrained to be 0.37!$ H^l '53. 
Interestingly, the region vv(grow) > w(geom), which should be generically favored by theories that slow struc- 
ture formation relative to general relativity, is quite restricted by data already. We find w(grow) < -0.80 at 2a. 
As an example, the best-fit flat Dvali-Gabadadze-Porrati (DGP) model approximated by our parametrization 
lies beyond the 3a contour for constraints from all the data sets. 



I. INTRODUCTION 

Observations of distant supernovae (SNe), galaxies, clus- 
ters of galaxies, and the cosmic microwave background 
(CMB) have shown that, surprisingly, the cosmic expansion 
is accelerating. This reveals that fundamentally new physics 
is missing from our understanding of the universe fl]. 

The cosmic acceleration may arise either from "dark en- 
ergy," a mysterious yet presently dominant component of the 
total energy density, or from "modified gravity," a modifi- 
cation of general relativity (GR) on large scales. The first 
case includes, for example, Einstein's cosmological constant 
or quintessence, a dynamical scalar field yfl. The second 
case includes modifications of four-dimensional GR due to 
the presence of extra dimensions, scalar-tensor theories, and 
others HHH. 

Current efforts focus, within the dark energy paradigm, on 
improving the constraints on the dark energy density £Ide, its 
equation of state (EOS) w = P/p and its time evolution dw/da 
(where a is the scale factor), by using observational data that 
bear on geometrical distances and the growth of structure. As 
first emphasized by flg] and subsequently discussed by many 
others [7], GR predicts a definite relation between geometrical 
distances and growth which is generically violated by modi- 
fied theories of gravity. To the extent current data (that are 
sensitive to different combinations of geometry and growth) 
yield consistent dark energy constraints, one can interpret this 
as a confirmation of the dark energy + GR framework. The 
simplest dark energy model, the cosmological constant, has 
passed this kind of consistency test so far |8|]. 

In this paper, we sharpen the consistency test. Our method 



goes by the name of "parameter splitting" as proposed by [9, 
llOll . Let us illustrate the technique using the cosmological 
constant (A)-cold dark matter (CDM) model. Instead of fit- 
ting the suite of observational data with a single cosmologi- 
cal constant density parameter ft\ (in addition to, of course, 
other non-dark energy parameters), we fit them with two pa- 
rameters i7A(geom) and f2A(grow): one determining the ge- 
ometrical distances, and the other controlling the growth of 
structure. The conventional approach is to assume the two 
parameters are equal. Here, they are allowed to vary sepa- 
rately. We employ the Markov chain Monte Carlo (MCMC) 
technique ifTTIl to derive the marginalized constraints on both 
parameters. If the ACDM model is correct, these two pa- 
rameters should agree within their uncertainties. This tech- 
nique of splitting a conventional parameter into two "meta- 
parameters" can of course be applied to any other parameter. 
In this paper, we will consider the splitting of both Oa and w. 

It is important to emphasize that parameter splitting checks 
for consistency not only between different data sets, but also 
for internal consistency within any single data set that is sen- 
sitive to both geometry and growth. In some sense, the con- 
ventional approach of obtaining constraints on, e.g., J7a sepa- 
rately from SNe, CMB, lensing and so on, and checking that 
they are consistent, is itself a simple form of parameter split- 
ting, i.e., splitting into SlA(SNe), SIa(CMB), fUOensing), 
etc. The parameter splitting that we employ here represents 
a more stringent, and theoretically better motivated, consis- 
tency test. It is also useful to note that there is a wide vari- 
ety of modified gravity theories. Our splitting of tt\ and w 
can be thought of as a crude way to parametrize the space of 
such theories. For instance, in the Dvali-Gabadadze-Porrati 
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(DGP) theory where gravity becomes weaker on large 
scales, structure growth is slowed and therefore one expects 
qualitatively w(grow) > w(geom) I6l [l2ll . 

We caution that should an inconsistency be discovered via 
parameter splitting, modified gravity is not the only possible 
interpretation. Systematic problems with the data, as well as 
complications in the dark energy model (such as a time vary- 
ing w or nontrivial dark energy clustering 1U3I0 . are also pos- 
sible. Additional parameters need to be introduced to check 
for the latter case. Parameter splitting can be applied to the 
new parameters as appropriate. 



II. GEOMETRY 

All geometrical distances in cosmology, such as the lumi- 
nosity or angular diameter distance, are related to the radial 
comoving distance 



X(z) = 



dz' 



(1) 



Jo H(z') 

setting the speed of light c = 1. The Hubble parameter H 
as a function of redshift z, i.e., the expansion history, can be 
parametrized as follows: 

H 2 (z) 



Hi 



= n m ( i + z) 3 + Q r ( i + z) 4 + n DE ( i + z) 



3(l+w) 



(2) 



where Ho = 100/z km s _1 Mpc _1 is the Hubble constant today. 
Throughout this paper, we assume that the universe is spatially 
flat, the dark energy has a constant EOS parameter w and all 
three species of neutrinos are massless. tt r is the radiation 
density today, in units of the critical density, including pho- 
tons and massless neutrinos; SIde is the present dark energy 
density, denoted as fl\ for the cosmological constant model 
(w = -1). Note that for a flat universe, the dimensionless mat- 
ter density Vl m can be replaced by 1 - Q r - Ude- We will use 
a superscript "(geom)" to denote the dark energy parameters 
appearing in the expressions of geometrical distances. 



III. GROWTH 

Inhomogeneities grow under gravitational instability ac- 
cording to the prevailing structure formation paradigm. The 
dynamics within the GR framework is described by a set of 
Boltzmann-Einstein equations well documented in the liter- 
ature 11411. In this paper, we use the publicly available code 
CAMB [ 15] to evolve these equations. For the purpose of illus- 
trating our method, and purely for this purpose, let us consider 
the special case of subhorizon matter fluctuations in the late 
universe. They evolve according to 5 m + 2HS m = AnGp m 5 m , 
where 8 m = 5p m /p m is the matter overdensity, p„, is the av- 
erage matter density, G is the Newton constant and the dots 
denote proper time derivatives. We ignore the dark energy 
perturbations here for simplicity. The growth equation can be 
rewritten as 



d S m 
din a 2 



d\nH 
dlna 



■ + 2 



d5„, 
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dlna 2a 3 // 2 
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where a = 1/(1 +z) is the scale factor. Therefore, the expan- 
sion history [Eq. (0] that determines geometrical distances 
also determines the growth of structure, in a way that is 
uniquely predicted by GR. 

It is not surprising that, in order to match existing data, vi- 
able theories of modified gravity often predict an expansion 
history (and therefore geometrical distances) that is similar to 
the one in Eq. (0. Such theories, however, generally predict a 
relation between expansion history and growth that is differ- 
ent from the one in Eq. (0. Given the wide variety of these 
theories, and in the absence of a particularly compelling can- 
didate [16], a crude way to test for such a possibility is to 
allow the dark energy parameters to take different values in 
the growth equation [Eq. (0] from their values in the expres- 
sion for distance [Eq. (Q])], i.e., parameter splitting. We use 
a superscript "(grow)" to denote the dark energy parameters 
characterizing the evolution of inhomogeneities. 

Note that one has some freedom in exactly how the param- 
eter splitting is performed. For instance, in Eq. (0, the dark 
energy parameters show up in two places: the second term on 
the left hand side of the equation (dlnH /dlna) and the term 
on the right hand side (fl m /H 2 ). One could choose to assign 
all of them to the "growth" category which is what we do, or 
one could assign some to the "geometry" category and the oth- 
ers to the "growth" category. Ultimately, there are many pos- 
sible consistency tests, and here we have chosen to perform 
one that is particularly simple to implement, i.e., assigning all 
dark energy parameters that enter the fluctuation equations to 
the "growth" category. It is worth noting that in a lot of mod- 
ified gravity theories, the equivalent of the Poisson's equation 
is often modified without modifying energy-momentum con- 
servation. In that case, one could argue assigning the term on 
the right hand side of Eq. (0 alone to the "growth" category 
might make more sense. We hope to investigate this in the 
future. 

The exact Boltzmann-Einstein equations for the evolution 
of structure, allowing for multiple components, photons, neu- 
trinos and so on, are more complicated than Eq. (0. The 
same parameter-splitting scheme can nevertheless be applied 
to the exact equations, which is what we do. This means, 
for example, the shape of the transfer function, such as the 
radiation-matter equality peak of the power spectrum, is de- 
termined by the growth parameters - recall that the transfer 
function is completely determined by the dynamics of fluc- 
tuation growth. The conversion of a feature, such as the 
radiation-matter equality length scale to an observed angle, 
on the other hand, involves the geometry parameters. 



IV. THE PARAMETER-SPLITTING TECHNIQUE 

To illustrate how the splitting of dark energy parameters 
into the "geometry" and "growth" categories is done in our 
analysis, we start with the weak lensing (WL) observables. 
There exists a natural division between the two categories for 
each term involved in the calculation @]. 

WL surveys measure the aperture mass statistic on different 
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angular scales 8: 



£d£P K (£)W 2 (£d), 



(4) 



where W is a window function with no dependence on cos- 
mology. P K {£) is the convergence power spectrum at the an- 
gular wavenumber £, given by 



£(z) 



-n 2 m H* I dz(i+zf 



d X (z) 



dz «gal(z') 



dz 

X(z')-x(z) ~ 
x(zO . 



X(z) 



(5) 



Here -Pi[^/x,z] is the matter power spectrum at wavenumber 
k = £/x and redshift z, n ga i is the normalized redshift distribu- 
tion of the background galaxies, and we have used Limber's 
approximation. We express everything in terms of the redshift 
Z, which is an observable of the surveys. 

Consider for instance the splitting of fl\ for the flat ACDM 
model. The three-dimensional matter power spectrum P$ and 
the mean matter density Q, m (= 1 - where the contribution 
of radiation is neglected at low redshifts) sitting outside the in- 
tegral both describe the foreground inhomogeneities through 
which photons travel. Therefore they go into the "growth" 
category and are calculated using 

fiferow), AU X ' S within the 

integral fall naturally in the "geometry" category. This in- 
cludes the x m the wavenumber £/x, which reflects the con- 
version between the observed angle and the physical length 
scale. These geometrical distances are all calculated using 
^(geom) ^ similar split can be applied to w in the context of 
the quintessence (Q)-CDM model. 

With the WL example in mind, we next consider the CMB. 
The temperature anisotropy power spectrum is given by 



C[ T = - I k z dkPy(k) 

7T 



e e (k, z = Q) 



*(*) 



(6) 



where ^(k) is the primordial metric perturbation (in confor- 
mal Newtonian gauge), Pyjk) oc k n ~ A is the power spectrum 
of ty, and &£ is given by 



e e (k,z = 0) = / dz' S T (k,z!)Mkx(z!)l 



(7) 



where ji is the spherical Bessel function and St denotes some 
source function. All the complicated dynamics is contained in 
St- Publicly available Boltzmann codes lfl5[[l7ll can be used 
to compute St, and therefore 9^, for any given primordial 
perturbation (Qt/ty is independent of "P; see UM ). 

We perform the geometry-growth split of Eq. (|7]i as fol- 
lows iflSil : St falls under the "growth" category and the rest 
(namely x(z') in the argument of jg) falls under the "geom- 
etry" category lfl9ll . The rationale for this particular way 
of splitting is most transparent when considering the Sachs- 
Wolfe term [ 20], where ST{k,z) is well approximated by Sd(z- 
Z*)[6o + \f r ](£,z«). Here #d(z— Z*) is the Dirac delta function 
with z* being the redshift of last scattering, and <do(k,z*) and 



^(fcjZ*) are the temperature monopole and metric perturba- 
tions at last scattering. Therefore, the Sachs- Wolfe term is 



0)~[Q + *](k,z*)Mkx*), 



(8) 



and our geometry-growth split is equivalent to using the 
growth parameters to compute [8o + ^](k, z*) and the geome- 
try parameters to compute x*, the distance to last scattering. 

It is straightforward to generalize the above splitting 
scheme to similar expressions describing the polarization 
spectrum. In the case of SNe, parameter splitting is trivial 
since SNe constrain only the geometry parameters. The split- 
ting for galaxy clustering is done as follows. As discussed 
earlier, the growth (as opposed to geometry) parameters de- 
termine the transfer function for the matter power spectrum. 
On the other hand, to measure the three-dimensional power 
spectrum of galaxies as a function of comoving spatial scale, 
one has to adopt a cosmological model in order to convert 
the observed redshifts and angular separations into comoving 
distances. This conversion is trivial for low-redshift surveys 
(involving only Ho) such as the Two-Degree Field Galaxy 
Redshift Survey (2dFGRS), but is nontrivial for moderate red- 
shift samples, such as the luminous red galaxies (LRGs) in 
the Sloan Digital Sky Survey (SDSS). For the LRGs, we fol- 
low luH] an d include a cosmology-dependent rescaling of the 
£-axes [22]. This rescaling is included in the "geometry" cat- 
egory. 



V. CURRENT OBSERVATIONS 

Below we list the four data sets used in our analysis. Many 
of these, though not all, are included in the CosmoMC pack- 
age dl. 



A. Cosmic Microwave Background 

We use (i) the recent Wilkinson Microwave Anisotropy 
Probe (WMAP) three-year data set and (ii) small scale 
CMB observational data including Arcminute Cosmology 
Bolometer Array Receiver (ACBAR) |0], Balloon Obser- 
vations Of Millimetric Extragalactic Radiation and Geo- 
physics (BOOMERanG) [25] and Cosmic Background Im- 
ager (CBI) (2(1 ■ We modify the Boltzmann code CAMB JH 
by splitting the dark energy parameters as described above. 
We assume adiabatic initial fluctuations, and neglect B-mode 
polarization and tensor modes. 



B. Supernovae 

We use the SNe data set for the Supernova Legacy Survey 
(SNLS) analysis described in H3]. 
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FIG. 1: Joint constraints on f^ eom> 



and Q. 



(grow) 



in a ACDM model 



(upper panel) and the normalized likelihood distribution of AQa = 
Qtgeom) _ £2<g rmv ) Qower panel). Here the equation of state parameters 
are fixed as vv <gcom) = w (gmw) = -1. The contours and curves show the 
68% confidence limits from the marginalized distributions. The thick 
gray lines show f^ com) = f^ low) . The data sets used are described in 
the text. Different contours and curves represent constraints from 
different combinations of the data sets. The smallest contour and 
the most narrow curve (black solid line) represent constraints from 
all the data. No significant difference is found and deviations are 



constrained to Af2/> 



-0.0044 



+0.0058+0.0108 
-0.0057-0.0119 



(68% and 95% C.L.). 
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FIG. 2: Variations of CMB temperature power spectra due to differ- 



ent changes of f^ com) and f2^ IOW ' (with all the other cosmological 



parameters fixed) as illustrated in the inset on the Sij 



vs. sZ; 



plane. The black solid curve corresponds to the black square symbol, 



which is our best-fit ACDM model with Q l f 



(gcom) 



(grow) 



0.744. 



The blue dashed curve corresponds to the blue circular symbol, 
which is obtained from the best-fit model by fixing f^ c ° ra) = f2^ row) 
and increasing both parameters by 0.03. The red dot-dashed curve 
corresponds to the red triangular symbol, which is obtained by fixing 
Qa and increasing f2^' ow) by 0.03 while decreasing f^ com) by 0.03. 



scales between 1 arcmin < 6 < 1 deg. To utilize the WL mea- 
surements on small scales, we take into account nonlinear ef- 
fects using (i) in the ACDM case, the nonlinear power spec- 
trum based on the halo model |32ll; or (ii) in the QCDM case, 
the mapping prescription in 113311 . 



C. Galaxy Clustering 



VI. ESTIMATING LIKELIHOODS 



We use data sets from (i) the Two-Degree Field Galaxy 
Redshift Survey (2dFGRS) which probes the galaxy dis- 
tribution at redshift z ~ 0.1 and the power spectrum on scales 
of 0.022/1 Mpc" 1 < k < 0. 1 8/i Mpc" 1 , and (ii) the luminous red 
galaxies in the Sloan Digital Sky Survey (SDSS) JH, which 
are at an effective redshift of z ~ 0.35 and cover scales be- 
tween 0.012/z Mpc" 1 < k < 0.20/1 Mpc" 1 . Redshift-space dis- 
tortions, galaxy biasing and nonlinear clustering 12911 are dealt 
with in ways described in lf2ll l28ll . 



D. Weak Gravitational Lensing 

Cosmic shear, due to weak lensing (WL) by large scale 
structures, has been detected by several groups [30]. The data 
set used in our analysis is from the 75 deg 2 Cerro Tololo Inter- 
American Observatory (CTIO) lensing survey 13111 . It covers 



We use the MCMC package CosmoMC ||23J to perform 
our likelihood analysis. CosmoMC uses CAMB [15] to cal- 
culate the temperature, polarization and matter power spec- 
tra. We modify both the CAMB and the MCMC por- 
tions to implement the parameter-splitting technique. In 
addition to the dark energy density and EOS parameters 
(n^ eom) ,f7*' ow, ,w (geom) ,w (grow) ), our cosmological parameter 
space includes the baryon density, the Hubble constant, the 
reionization optical depth, the scalar spectral index and am- 
plitude of the primordial power spectrum: (ilbh 2 ,h,T,n s ,A s ). 
When w (grow) 1 , sound speed of the dark energy is set as 1 in 
CAMB [ 15]. For simplicity, we assume a flat universe for both 
geometry and growth parameters. The Monte Carlo chains 
are generated by the Metropolis-Hastings algorithm ll34ll . We 
adopt Gaussian priors of Sl^/i 2 = 0.022 ± 0.002 from Big Bang 
nucleosynthesis (BBN) 01 and H = 72 ± 8 km s^Mpc" 1 
from the Hubble Space Telescope (HST) key project [36]. 
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VII. RESULTS 

Applying our consistency test to the ACDM model, where 
the EOS parameters are fixed as w (geom) = w (grow) = -1, the 
upper panel in Fig. Q] shows the marginalized constraints on 
plane. The confidence contours follow 



the n { l mv,) vs. n ( l eom) 

roughly, but not exactly, the f2' 



(geom) 
A 



n 



(grow) 



line. The inter- 

(geom) 



esting quantity in this case is the difference ASIa = f2 A 
f^ row) , wnose normalized probability distribution is shown in 
the lower panel of Fig.Q] When all data are utilized, we find 
the marginalized constraint Afl A = -0 .0044^ 5 5 ^° \% (68% 
and 95% C.L. respectively). Figure Q] also shows that CMB 
anisotropics, when combined either with galaxy clustering or 
SNe, deliver most of the overall constraining power, i.e., hav- 
ing the narrowest likelihood distributions. 

We also find the marginalized constraint on the aver- 
age Tt A = (Q ( l eom) + nf r0VJ) )/2 using all data sets: Ti A = 
0.744^;o[g+oo3i. The constraint on the difference is almost 
three times better than the constraint on the average. The 
CMB contour in Fig. Q] even without the addition of other 
data, already exhibits this trend. Let us therefore focus on 
understanding this phenomenon in the context of CMB. 

As illustrated in Fig. [2] increasing both f^ 60 ™' and f^ 1 ™^ 
by the same amount (with all the other cosmological pa- 
rameters fixed) produces only a small shift of the predicted 
Ci (blue dashed curve). However, moving in the orthogo- 
nal direction, i.e., increasing f2^ low) while decreasing 51^ eom) , 
creates a much larger shift (red dot-dashed curve). It ap- 
pears partial cancellations occur between the shift in the 
distance to last scattering (a geometrical quantity) and the 
shift in the sound horizon (which controls fluctuation growth) 
when one changes both 51^ eom) and ^! A !1C,W, by the same 
small amount, creating a roughly degenerate direction along 
fi (geom) = ^(giow) Conversely) the effects of the two different 

shifts roughly add when one changes il A ;ec,m ' ) and f^ 1 ™^ in 
opposite directions, making Ail A highly constrained. 

One could argue that in theories of modified gravity con- 
structed to explain the late time cosmic acceleration, the 
growth of fluctuations should only deviate from GR at late 
times. A better approximation of such theories is perhaps 
to split the EOS parameter w. We therefore next apply our 
consistency test to the more general QCDM model. The 
EOS parameters, w (grow) and w (geom) , are assumed constant, 
but are allowed to vary independently. In this test, we as- 
sume riDE° m) = i!pg OW \ The upper panel in Fig. [3] shows the 

(geom) 



plane [37]. 



marginalized constraints in the w (grow) vs. w 1 
We again find that the difference Aw = w (geom:) -w (grow ' 1 is con- 
sistent with zero; deviations are constrained by combining all 
data to Aw = 0.37^q3^q 53 (lower panel in Fig.0. The average 
is constrained to be w = (w (geom) + w (grow) ) /2 = -1.1 3^;^g;g. 

Figure |3 shows a long tail towards large negative values of 
w (grow)^ wn j cn can b e understood as follows. Density perturba- 
tions can grow significantly only during the matter-dominated 
epoch, and as w (glow) becomes more negative, this epoch is 
longer (i.e., dark energy domination occurs more recently). 
The extension of the likelihood contours in the large negative 
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FIG. 3: Joint constraints on w <gcom) and w (grow) in a QCDM model 
(upper panel) and the normalized likelihood distribution of Aw = 
w (geom) _ w (grow) (j ower p ane [y Here the energy density parameters 



are fixed as Q 



(geom) 



.-.(grow) 
DE • 



The contours and curves show the 



; confidence limits from the marginalized distributions. The thick 
gray lines show w <geom) = w (glow) . The data sets used are described in 
the text. Different contours and curves represent constraints from 
different combinations of the data sets (see legend in Fig. [TJ. The 
smallest contour and the most narrow curve (black solid line) repre- 
sent constraints from all the data. No significant difference is found 
and deviations are constrained to Aw = 0.37^f 6 ^f 3 (68% and 95% 
C.L.). The star-shaped symbol corresponds to the effective w ge 
and w (grow> , which approximately match the expansion history and 
the growth history, respectively, of a flat DGP model with our best- 
fit Q,„. 



direction of w ( - gmw - ) reflects the fact that a very recent dark en- 
ergy domination is actually acceptable as far as the growth 
of structure is concerned. This does not imply the data is 
consistent with the absence of dark energy, however. On the 
contrary, the data prefer a low il„, which for a flat universe 
implies the presence of SIde- It is interesting to note that qual- 
itatively, the DGP theory prefers w (grow) > w (geom) a 
region that is quite restricted by data already. In fact, we find 
that a DGP model with our best-fit Q m , represented effectively 
by the star-shaped symbol in Fig. [3] lies beyond the 3ct con- 
tour for constraints from all the data sets; varying fl m in the 
DGP model within its 3er limits has little effect on the position 
of the point. We also find the upper limits of w (grow) < -0.97 
at lcr and w (grow) < -0.80 at 2a Q. 
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VIII. DISCUSSIONS 

Our study reveals no evidence of a discrepancy between the 
two split meta-parameters. The difference is consistent with 
zero at the la level for the ACDM model and 2a level for 
the QCDM model. We find tight constraints from the exist- 
ing data sets, especially on the difference between SI a derived 
from growth and £7 a derived from geometry (better than 1%). 
In other words, the cosmological constant model fits current 
data very well. Current data do not appear to demand mod- 
ified gravity theories. Parameter splitting can be thought of 
as a crude way to parametrize the space of these theories. As 
such, our constraints can be viewed as putting restrictions on 
modified gravity theories, but the precise constraints on any 
particular theory must be worked out on a case by case basis. 
The kind of constraints we obtain here are likely to signif- 
icantly improve in the future, as the cosmological data im- 
prove in quality and quantity. The power of future surveys is 
demonstrated by a calculation that a Large Synoptic Survey 
Telescope (LSST)-like survey could constrain Aw to 0.04, 
using shear tomography alone, an order of magnitude better 



than current constraint from all data sets l9l [39tl . 
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